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We calculate the time-evolution of a discrete-time fragmentation process in which clusters of particles break 
up and reassemble and move stochastically with size-dependent rates. In the continuous-time limit the process 
turns into the totally asymmetric simple exclusion process (only pieces of size 1 break off a given cluster). We 
express the exact solution of master equation for the process in terms of a determinant which can be derived 
using the Bethe ansatz. From this determinant we compute the distribution of the current across an arbitrary 
bond which after appropriate scaling is given by the distribution of the largest eigenvalue of the Gaussian unitary 
ensemble of random matrices. This result confirms universality of the scaling form of the current distribution 
in the KPZ universality class and suggests that there is a link between integrable particle systems and random 
matrix ensembles. 

S 1 ' 

. I. INTRODUCTION 

Asymmetric exclusion processes are paradigmatic models for systems far from equilibrium, both for their wide range of 
applications and the availability of exact results |ll|2|]. Despite their greatly reduced complexity they capture various generic 
nonequilibrium phenomena such as the occurrence of shocks and particle condensation ylHo!- A major breakthrough in the 
exact calculation of universal properties came with the realization that various ensembles of random matrices occur in the study 
of current fluctuations and related quantities in the totally asymmetric simple exclusion process (TASEP) |6, 7] and also in 
the polynuclear growth model 

MM El. 

Here we use random matrix theory combined with the Bethe ansatz to study a 
process introduced by Priezzhev 1 12] which describes the stochastic fragmentation and reassembly of diffusing particle clusters 
in discrete time. 

In the totally asymmetric fragmentation process fl2Tl one considers a one-dimensional lattice where each lattice point is 
occupied by at most one particle. A string of n consecutive particles (and bounded by empty sites) is considered a cluster of size 
n. The stochastic time evolution occurs in discrete time steps. From each cluster a piece of size n' < n may break off and move 
to the right by one lattice unit with probability (1 - p)p"' . The whole cluster may hop with probability p" . According to this 
definition the transition rules are the following: 



o 
o 



> 



OAA . . . AAA -> AA . . . AA OA with probability p( 1 - p) 



n n-l 

^ ; AA.^AAA 0^0 AA^A OAA with probability p 2 ( 1 - p) 

: " " 2 

o . 

OAA... AAA -> 0A0 A . . . AAA with probability p n ~ \l-p) 

£^ n n-l 

^ ! OAA . . . AAA -> 00 AA .. . AAA with probability p" 

Ch ! n n 

o . 

Notice that two clusters can merge through hopping if they are separated by only one vacant site. Other fragmentation processes 
have been studied recently 1 13l ll4l.ll5l.llal . 
• i-H , In the noiseless limit p — 1 neither fragmentation nor recombination occurs and all clusters move ballistically with probability 
1. In the limit p — > with an appropriate rescaling of time (t — > oo with t' - pt fixed) only pieces of size 1, i.e., single particles 
may break off. This process is equivalent to the usual TASEP. In discrete time the fragmentation process may be interpreted 
as an asymmetric exclusion process with long-range hopping. This is forbidden in the usual discrete-time TASEP with parallel 
update studied in Ref. |6]. Moreover, the fragmentation process has no particle-hole symmetry and the stationary distribution is 
a product measure with constant density p and stationary current 
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The absence of correlations implies that fragmentation and recombination balance each other such that no condensation into 
macroscopic clusters occurs. Indeed, our interest is not in the stationary state of the system but how it evolves from a fully 
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"phase-separated" initial state where the left half of the infinite lattice (all sites k < 0) is occupied while the right half (all sites 
k > 0) is vacant. 

Under Eulerian scaling (lattice constant a and time step t tend to zero with the ratio a/r kept constant) one expects the 
coarse-grained density p(x, t) to be governed by the hydrodynamic conservation law 

dtp + d x j(p) = (2) 

with the macroscopic current ©■ Our interest is in the microscopic fluctuations of the current which after appropriate rescaling 
is expected to be given by an universal scaling function. In particular, one obtains universal corrections to below the Euler 
scale. 

In Sec. 2 we first consider the continuous-time limit of the model, i.e., the TASEP. We present a new derivation of the result 
Eq. (1.18) of Ref. |6] for the exact distribution of the time-integrated current. We do not use the combinatorial arguments 
employed by Johansson, but show how the same expression can be obtained logically independently from the Bethe ansatz 
solution via the determinant expression of Ref. fnll . In Sec. 3 we extend this approach to calculate the current distribution for 
the fragmentation process with arbitrary fragmentation parameter p. In Sec. 4 we present some conclusions. Some properties of 
the functions we use in the calculation are given in the appendices. 



II. CONTINUOUS TIME TASEP 



A. Known results 



1. Current fluctuations 



Consider the TASEP on an infinite chain in continuous time where particles hop to the right with rate 1, provided the target 
site is vacant. At t — the left half of the system (from site -oo to site 0) is occupied while the right half is empty. We focus 
on the probability P(M,N, t) that the Mh particle (which was on the (1 - AOth site of the infinite cluster at t — 0) hops at least 
M times up to time t. Using combinatorial arguments involving the longest increasing subsequences of random permutations 
Johansson has proved (Eq. (1.18) in Ref. [6]) 

P(M, tf, t ) = _L f f[ ( Xi - xj? f[ xf V ',/\v (3) 

z m,n Jump i<i<j< N j=l 

for M > N. To have proper normalization the partition sum has to be 

Z'mn= f I! (Xi-Xj) 2 f]xf- N e^d N x. (4) 

The expression (|3j is equal to the probability that the largest eigenvalue of a random matrix AA* is < t where AisaNxM matrix 
of complex Gaussian random variables with mean zero and variance 

1/2 ill- 
Let J(x, t) be the number of particles that have crossed the lattice bond (x, x + 1) up to time t, i.e., the time-integrated current. 
By construction one has for the probability that J(x, f) > m 

Prob [J(x, t) > m] - P(m + x + 1, m + 1 , f). (5) 



2. Bethe ansatz 



Consider the TASEP on the infinite chain with a finite (AO number of particles located initially at sites A^ = Huh,'" , h) 
(l\ < h < ■■■ < In)- It is proved in 1 1 711 that the probability of having these particles on sites = {k\, k-i, ■ ■• , k^} (k\ < < 
• • ■ < kfj) at time t is given by the determinant 



Q{A N ,B N ;t) 

F~N-i(kN -h',i) F N -2(kN - h\ i) 
For the definition and some properties of the F p (n; t) functions see appendixlAl 



Foih-lut) 
F!(k 2 -lut) 



F-iih 
F (k 2 - 



-h;t) 
hit) 



F-N+i(ki ■ 
F-N+i(ki ■ 



In, t) 
In, t) 



Fo(k N - In', t) 



(6) 
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The determinant has been obtained from a coordinate Bethe ansatz for the conditional probability Q(A^, B^; t). Given the de- 
terminant, the proof that it is the solution of the master equation for the TASEP follows from standard relations for determinants, 
for details see 11711 . In the next subsection we show how (|3j can be derived directly from l|6} without reference to combinatorial 
properties of the process. 



B. Calculation 



The dynamics of the rightmost N particles in the TASEP are independent of all particles to their left. Therefore P(M, N, t) of 
sec. Ill A fl can be expressed via Q(A^, B^; t) of sec. Ill A 21 



P(M. N.n = V Q({-N + 1, -JV + 2, • • • , 0}, {Jki, Jfe 2 , • • • , k N ); t). 

M-N<k\ </(2<—<kN 



(7) 



Inserting Q and (|6} one gets 



1 r N 

w jbwf i<T</<iv y 



z 



M-N<k\ <k 2 <--<kf, 



F (h+N-l;f) F- X (ki+N-2;i) 
Fdki + N-ht) F (k 2 +N-2;t) 

Fat-i (k N + N - 1 ; f) F N - 2 (k N +N-2; t) 



F-u+iihU) 
F- N+2 (k 2 ; t) 

F (k N ;t) 



(8) 



In what follows we show by determinant manipulations and using properties of the functions F^ that this equality holds. 

Our starting point is the rhs. of (|8j. The summation over (ki,k 2 , ■ ■ -k^) can be done in N steps for which we choose the 
following sequence: 



k N -i 



(9) 



z =z z - z z 

M-N<k\<k 2 <-<k N k N =M k N _ l= M-l k 2 =M-N+2 k^M-N+l 

After summation over k\ (for which we use (IA5» the first row of the matrix becomes 

Fi(M; t) - Fi(k 2 +N- F (M - 1; t) - F (k 2 + N - 2; f) ••• F_ N+2 (M - N + 1; t) - F_ N+2 (k 2 ; t) (10) 

which reduces to 

F x (M;t) Fo(M-l;t) ■■■ F. N+2 (M - N + 1; t) (11) 

after adding the second row to it. The same method can be used up to the sum over kff-i. For the last sum we use dA7> and 
finally we get 



Fi(M; t) 
F 2 (M+l;t) 



F (M-l;t) 
Fi(M; t) 



F N (M + N-l;t) F N _i(M + N-2;t) 



F- N+2 (M-N + l;t) 
F- N+3 (M~N + 2;i) 

F\(M;t) 



(12) 



for the rhs. of i). 

It turns out to be useful to represent all the F functions in \\2\ by an integral using (IA4> . Since for M > N all are zero at t — 
the determinant (II 21 can be written as 

Jf F (M - l;r)dr jf F^(M -2;r)dr ••■ jf F- N+1 (M - N;r)dT 

^F^M-^dT £F (M-l; T )dT ■■• fiF- N+2 (M-N + l;z)dT 

jf F N -\(M + N - 2; T)dT jf F N . 2 (M + N-3; T)dr ■■■ jf F (M - 1 ; t)</t 
In the second row we perform a partial integration after which the 2th element becomes 



(13) 



tF 2 -i(M -i+l 



Jo 



tFi_,-(M - i; r)dr 



(14) 
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Note that the constant part is t times the corresponding element of the first row so after substracting t times the first row the 
second row becomes 



tF (M- 1;t)c/t 



f 

Jo 



tF-i{M -2;r)dT 



f 

Jo 



tF-n+\(M - N; r)dr. 



(15) 



In the third row we perform a double partial integration then we add f 2 /2 times the first row and substract t times the (original) 
second row. Repeating the same procedure for all of the rows we get 



J£f (J#-1;t)</t 
~~ Jq tF o(M - 1 ; T)dT 

\lT 2 F Q {M-\-T)dT 
±£T 3 F (M-l;T)dT 



j£F_i(Af-2; T)dT 
-£tF_ 1 (M-2;t^t 
\] Q T 2 F^{M-2-r)dT 
-hJ^ F -iiM-2;T)dT 



Jq F-n+i (M — N; r)dr 
-CTF_ N+1 (M-N;r)dT 
\] Q T 2 F- N+x {M-N-T)dT 
-jifoT 3 F_ N+l (M-N;T)dT 



(N- 



^{V-'^M-l;^ _^_jrV- 1 F_ 1 (M-2;r)dT •■■ ^ j^F^CM - N;r)dj 



^ '! J[0/]« 



F (M- l;n) F_i(Af - 2; n) 

F (M-1;t 2 ) F_!(M-2;t 2 ) 

F (M-1;t 3 ) F_!(M-2;r 3 ) 

F (M-l;r w ) F_i(M-2;r^) 



F_ N+l (M-N; Tl ) 
F. N+1 (M - N;r 2 ) 
F_ n+1 (M-N;t 3 ) 



F- N+1 (M-N;t n ) 



d N T. (16) 



Using JA6> we can rewrite F_i(M - 2;t,) as Fq(M - 2;t,) - Fq(M - l;r,). So adding the first column to the second one 
the zth element of the latter becomes Fq(M - 2; r,). Similar transformations can be done with the other columns as well: since 
according to ( IA8I 



F p {n;t)=Y J {-\)4 P )F (n + m;t) for p < 0, 

m=0 ' m ' 



(17) 



one can add suitable linear combinations of the first I - 1 columns to the Ith one so that the zth element of the fth column becomes 



F (M-/;r i ) = e- 



(M - /)! 



(18) 



After these transformations one gets for (fTBl : 

JV-l , AT , r N 

(->) [ «nsn^Ln(^-")^-<-' 



T A<-1 ^JV-2 



-AM tJV-2 
'AT 'A 1 



T 1 T° 

J J 

T 2 T 2 



1 

'A< 'a? 



d t. 



(19) 



Since the integration is symmetric in the r, while the determinant is antisymmetric we may replace the product of the r, by 
antisymmetric combination 



1 

AM 



T? T 1 



i ' i 

"-2 _JV— 1 



'A 1 1 N 



-N-2 _JV-1 
'A 1 'A? 



(-!)[« 



AT! 



r , , N-2 
' 1 ' 1 

T N - 2 



N-\ N- 
l N 'N 



T 1 T° 
J 

T 2 T 2 



1 
1 AT N 



(20) 



This is the Vandermonde determinant 



T N-1 T N-2 

Jf-i J/-2 

T 2 T 2 



T N-1 N- 

'n 'n 



1 

T 2 T 2 



1 
N l N 



n ( t *- t a 

1 </< j<N 



(21) 
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and leads to 



/=1 X V ' 1 ^ 1U ' fJ i=l l<i<j<N 

This is in agreement with the lhs. of (|8j, moreover we get the partition function as a "by-product": 

N 

z M,iv=n /!(M "° ! < 23) 



C. Asymptotic form of the density profile 

For the TASEP the current Q reduces to j = p(l - p) and the density profile of the step-function initial state evolves on Euler 
scale according to 

p(x,f) = i(l-v) (24) 

where v = x/t and -1 < v < 1 lfl9ll . Outside this range the density keeps it initial value. 

To describe the density profile below Euler scale we use the asymptotic form of the result l|5} for the current distribution. It 
can be shown that that the formula ( 1221 can be written as a Fredholm determinant with the Meixner kernel which in a proper 
limit reduces to the Airy kernel (for details see |6]). This implies that in this limit the asymptotic form of P(M,N, t) is given by 
the Tracy-Widom distribution 1 20] of the Gaussian unitary ensemble (Fq UE ), viz. 

lim P([yN],N, a>(y)N + cr(y)N l/3 s) = F Qm (s) (25) 

N—>co 

with 

io{y) = (1 + Vr) 2 and cr(y) = y~ 1/6 (1 + Vr) 4/3 ■ (26) 

This type of scaling is characteristic for one-dimensional surface-growth models of the KPZ universality class (|8|). 

For our initial state the integrated current J(x, t) gives the number of particles being at sites k > x at time t. It can be shown 
using theorem 1.6 of 1 6] that 



lim Prob 

t— >oo 



J([vt],t)<Ul-vf + 2-^(l- V f 3 t^ S 



1 - F CI]E (- S ) (27) 



for < v < 1 . Because of particle-hole symmetry one has a similar expression for - 1 < v < 0. 
We can also consider the corresponding surface growth model where h(x, f) is defined as 

h(x, t) = \x\ + 2J(x, t) (28) 

The asymptotic mean shape of h([vt], t)/t = (1 +v 2 )/2 on Euler scale follows directly from 11241 . The deviations can be calculated 
from i21\ : 



lim Prob 

t— >oo 



^(l + v 2 )-/ ; ([v f ], f) <2- 1 / 3 (l-v 2 ) 2/3 f 1 / 3 , 



fW*) (29) 



for < v < 1 . This implies that 

U\+v 2 )-h{[vt],t) 



= -c, (30) 



2-1/3(1 _ v 2)2/3 f l/3 

where h(x, t) is the mean height at position x and time t, and 

dssF' Gm (s) = -1.77109 (31) 

CO 

is the mean of the distribution Fgue- Since the density p(x, t) is |(1 - h(x, t) + h(x - 1, f)) expression ( l3*0l allows us to calculate 
the correction to the density profile J24i 

p( vt ,t) = —-+c— v(l-v 2 ) r 2 ' 3 . (32) 

The KPZ exponent 2/3 appearing here is universal. Note that the coefficient v 1 1 - v ) diverges as v — > + 1 . In these singular 
points we expect af 1 ' 2 correction. 
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III. THE DISCRETE-TIME FRAGMENTATION PROCESS 



A. Complete solution of the master-equation by Bethe ansatz 

The solution is very similar to the one of the continuous time TASEP. As in Ref. one constructs the Bethe solution for 
the master equation for the conditional probabilities and recasts the expression in terms of a determinant. Then one proves by 
elementary matrix manipulations that the analogue of equation ijfji is the following (for details see [ 12]): 



Doih-lrj) D- X (k x - l 2 ;t) ■■■ D- N +i(ki-l N ;i) 

Di(kz-h;t) D () (k 2 - h;t) ■■■ D- N+ z(k 2 - l N ;t) 

Q(A N ,B N ;t) = 

D N -\{k N -l x ;t) D N _ 2 (kN~h;t) ••• D (k N -l N ;f) 
For the definition and the main properties of the D functions see appendixlBl 



(33) 



B. Calculation of P(M, N, t) 

We follow the strategy of the previous section, but some care needs to be taken when manipulating sums rather than the time 
integrals. Here P(M,N,t) is again defined by ([7) but now with the Q of (1331 . The summation gives the same result as for the 
continuous time version since the corresponding properties of the D functions are the same as those of the F functions: 



P(M,N,t) = 



Di(M; t) 
D 2 (M+l;t) 



D () (M-l;t) 
D X {M; t) 



D N (M + N-l;t) D N ^(M + N-2;t) ■■■ 
We can represent each element by a sum over the discrete time variable if M > N: 



D- N+2 (M- N+ l;r) 
D- N+3 (M-N + 2;t) 

Di(M;f) 



P(M,N,f) = p h 



Z£m-iA)(M-1;0 



2;r = 1 M _ 2 £»-i(M-2;f) 



Z>= M Di(M;t') 



2TJm-iA>(M-1;0 



2Xm+jv-2 d n-i(M + N - 2; f) Z'^ M+N -3 d n-2(M + N-3;f) ■■■ 
Instead of the partial integration we perform partial summation here and use the identity 

t-i t-i 



2T=m-iA)(M-1;0 



^ bAa.f + i - a,,) = ~2_j at+iibr+i - b t >) + a t b, - a h b h . 



i'=ii 



for the second row with a t = t — 1: 

t-i 



^Dy(X,t') = -pY^t'Dy^X - l,t') + (t- l)Dy(X,t)-(X- l)Dy(X,X) 

f=x t'=x 

t-i 

= -P^] t'D v ^{x-\,t') + {t-\)D v {x,t). 

f=x-i 

For the third row we go further and apply (I36> again with a, = t(t - l)/2: 

t-i 

- P ^ t'Dy-\{X -\,t') + (t- l)Dy(X, t) = 



f=x-i 



t'(t' + l) , f(f-l) (x-l)(x-2) 

2 D v - 2 (x - 2, t') - p^-t-Dy-tix -1,0+ P- \ -Dy-dx -l.x-l) + (t- l)D v (x, t) = 



2 V< t'(f + \) 



Z 

f=x-2 



(34) 



(35) 



(36) 



(37) 



, Dy^(x-2,t')-p t -^^-D v ^(x-l,t) + (t-l)Dy(x,t). (38) 
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In the fourth row we perform the partial summation once more with a, — (t - l)t(t + l)/3! and similarly for all the rows. Finally 
after adding suitable linear combination of the first n — 1 rows to the n-th one we get 



P(M,N,t) = 

P 2 (-1)1 2} 



0!1!---(W- 1)! 



ZjrVi A)(M - 1; f) Z>'m-2 D-\(M - 2; f) 
Z£a#-i t>D o(M - 1; f) t'D-i(M - 2; f) 



Z't~=M-i f N ~ 1 D (M -l;f) Z>=m-2 t' N - [ D^(M - 2; f) 
We can set all the lower limits of the sums to and write 



Z'-M-N D l-N(M-N;t') 

Z'^M-Nt'Di-NiM-N; f) 

Zr=M- N t' N ~ lD i-N( M -N;t') 



P(M,N,t) 



p 2 (-l)LiJ 



0\1\---(N- 1) 



Z AA^, ■ ■ ■ 



h,ti,— ,tN=0 



D (M-l;h) D-i(M-2;n) ■■■ D^ N (M - N;h) 
D (M-l;t 2 ) D-i(M-2;t 2 ) ■■■ D X - N (M - N;h) 



D (M-l;t N ) D^(M-2;t N ) ••■ D^ N (M - N;t N ) 
Adding suitable linear combination of the first n—1 columns to the n-th one we get (by using jBll^ ) 



P(M,N,t) 



P 2 (-l) L2j 

0!1!---(N-1)! 



r-i 



/ j 'l'2 ( 3 l N 



«i,«a,-.«v=o 



D (M-l;h) D (M-2;h) ■■■ D (M-N;h) 
D (M-l;t 2 ) D (M-2;t 2 ) ■■■ D (M-N;t 2 ) 

D (M-l;t N ) D Q (M-2;t N ) ■■• D (M-N;t N ) 



Using the explicit form of Do JB8I one arrives at 
P(M,N,t) = 



p MN (l _ pr MN + ^l ( _ l) [i] 



t-\ 



0!1! •■■(#- 1)\(M- 1)\(M - 2)\ ■ ■ ■ (M - N)\ 

fi(fi - l)---(h -M + 2) h(h - l)---(h -M + 3) 
t 2 (t 2 - 1) ■ ■ -(h - M + 2) t 2 (t 2 - l)---(t 2 - M + 3) 



2 (l-p)^t\t\^..f N - 



h(h - l)---(h-M + N+l) 
t 2 (t 2 -l)---(t 2 -M + N+l) 



t N (t N - 1) ■ ■ ■ (t N - M + 2) t N (t N - 1) ■ ■ ■ (t N - M + 3) • • • t N (t N - 1) ■ ■ ■ (t N - M + N + 1) 



(39) 



(40) 



(41) 



p m (l-p)- 



u%iU-mM-jy. 



2 (1 - pfi* 'n\t\t\ ■ ■ ■ f N ~ l Y\ ~ W) -2)---(tj-M + N+l) 

h, h, -,fjv=0 j=\ 



fN-l fN-2 

'l H 

JV-1 t N-2 

l 2 l 2 



JV-1 fN- 
l N N 



p MN (l - p) 



-MN+m±<l t-1 N 

N , (M _T7- £ n^-lX0-2)-(^-M + iV + lXl-^ 

llj=l J-\. M .)>■ t u H,-,t N =0 j=l 



JV-1 fN-2 

'l 'l 

fN-l fN-2 

l 2 l 2 



JV-1 JV-2 



(42) 



Finally we get 



P(M,N,t) = 



p MN d-p) 



(M-N-l 



s n n*-» 



^ jfc=0 



(l-p)° 



(43) 



One can see that in the limit t = t/p, p —> 0,t = const. (I43t is equivalent to i22\ as expected. 
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-2 -1 012345 -2 -1 012345 

lattice site lattice site 



FIG. 1 : The first figure shows a possible history of three particles from time to time /. We adjust probability p to the diagonal lines (thick 
solid lines) which correspond to particle hopping. The weight of the thick dashed vertical lines is 1 - p while the thin solid vertical lines have 
weight 1. The second figure shows the corresponding path configuration of the DTASEP. The path of the second (third . . . ) particle is shifted 
by one (two . . . ) time steps. 

Expression J43i can be written as 

P(M,N,t)=Z{M,Ny l Yj II \Im-N Y 1 -*)'' n^-^ 2 - (44) 

t u h,- ,tf/=0 j=\ " ' ' i<j 

In the latter form P(M,N, t) is similar to the same quantity of the usual discrete time ASEP (DTASEP) (see Proposition 1.3 of 
01). Namely: 

P{M, N, t) = P D tasep(M, Nj + N-l) (45) 

This correspondence can be seen in a graphical representation of the dynamics (Fig. 1). The statistical weight coming from the 
path of the first particle is the same in both processes since its dynamics is identical in the two cases. The weight of the path of 
the Mh particle from to t is the same as that of the normal DTASEP from to t + N - 1 which explains the relation ( 145 \ (see 
also figure lilTBl . However, this correspondence is valid only in the case of this special initial condition. 

From the geometric interpretation it follows that ( 145 \ is valid for all values of M and N, although we derived this (on the 
level of formulas) only for M > N (note that Proposition 1.3 of 0] as well as d44b is valid only for M > N). The DTASEP has 
particle-hole symmetry which implies 

f D TASEp(M, N, f) = PdtasepW M, t), (46) 

from this we get 

P(M, N, t) = P D tasep(M, N,t + N-l) = Pdtasep (N, M, t + N - 1) = P{N, M,t + N- M). (47) 



C. Asymptotic form of the distribution function 

Knowing the derivation of the asymptotic form of Pdtasep(A^, N, t) (see section 3 of |6J]) it is rather easy to obtain similar 
results for this process. We repeat the result of Johansson for the discrete-time TASEP with parallel update with our notations: 



lim Pdtasep ( ly N], N, N aj DT ASEp(y, p) + N in o-(y,p)s) = F G ue(s) 



(48) 



9 



with 



(1 + V(wF) 2 

w D tasep(7> p) = + y, (49) 



(1 _ p)l/6 r -l/6 2/3 2/3 

o-(r )jP ) = - — £ — (Vr+Vw) (i+ Vd-p)r) • 



p 

This formula is derived for y > 1 in |6] but it is easy to show that it is valid also for < y < 1. Using ( I46l > one gets 



(50) 



lim p(yN,N,Naj DTASEP (l/y,p)y + N 1/3 cr(l/y,p)y' /3 s) = F Gm (s), (51) 

which is identical to d48l since wdtasep(1/t>p)t = wdtasep(t,p) and cr(l/y, p)y 1/3 = cr(y,p). 

To obtain the corresponding asymptotic result for the fragmentation process one goes through the same steps as Ref. 0] (or 
simply use (I45» and (I48> . One finds 

lim p([yN],N,No)(y,p) + N l/3 cr(y,p)s) = F cl]E (s) (52) 

with 

co(y,p)= + r -l = . (53) 

This result is valid for any y > 0. The corresponding identity is w(l/y, p)y — I + y - ai(y,p). We remark that amplitude cr of 
the deviation is the same for both models. 

IV. CONCLUSIONS 

We have shown that the distribution of the time-integrated current for the TASEP as well as for the totally asymmetric frag- 
mentation process can be obtained using Bethe ansatz. To show this one uses a determinant representation of the Bethe wave 
function which solves the master equation for a finite number of particles. After appropriate scaling the current distribution is 
given by the distribution of the largest eigenvalue of a random matrix ensemble. This observation may lead to better under- 
standing of the relation between the random matrix theory and the Bethe ansatz and suggests that also other integrable hopping 
processes may be treated in a similar fashion. The main task that remains is the derivation of suitable determinant representations 
for conditional probabilities for such processes. In the scaling limit of the fragmentation process the distribution of the time- 
integrated current around its mean converges to the same Tracy-Widom distribution for the Gaussian unitary ensemble found 
previously for the TASEP, thus confirming universality of this quantity. 
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Note added: After completion of this work we learned that the result of Sec. JII Bl has been found also by Nagao and 
Sasamoto, cond-mat/0405321 

APPENDIX A: F FUNCTIONS 

Definition: 

= ^-.& e-V-^'a-zT^dz. (A2) 
2m Jf z | = i_o 
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The integral is taken in a way that the pole at z = 1 should not be taken into account. 
Differentiation, integration: 

jF p {n;t) = F p ^{n-\\t) 



Summation: 



For p < and n > 0: 



f 2 F p (n;t) = 
Jt, 



F p+l (n + l,f 2 )- F p+ i(n+ l,h) 



^ F p (n;t) = F p+ i(ni;t) - F p+ i(n 2 + l;f) 
F p (n; f) = F p+l (n; t) - F p+l (n + l;f) 

CO 

^ F p (n;t) = Fp+iinut) 



F p (n;t)=e- t J](-iy 



m=0 



to j (m + n)\ 



(A3) 
(A4) 

(A5) 
(A6) 
(A7) 

(A8) 



F p (n; f) can be written as a (finite or infinite) sum also in other regions of the («, p) parameter space (see fl71) but those formulas 
are not used here. 

For n > (and any p): 



FJn; 0) = 5 nX 



(A9) 



Definition: 



APPENDIX B: D FUNCTIONS 



D q (n, t )=±.f*dk(l-p + pe- >k ) (l - e *+'°>) 
2m Xi = , n V Z I 



-P Jkn 



|z|=l-0 



Discrete time derivative and summation for t : 



Summation for n: 



D q (n, t + 1) - D q (n, t) = pD q - { (n -1,0 

^D 9 (n,f) = -{p q+ i(n +l,t 2 + l)-D q + x (n + l,h)) 

'='i p 



D q (n;t) = D q+1 (nut) - D q+1 (n 2 + l;f) 

n=n\ 

D q (n; t) = D q+l (n; t) - D q+l (n + 1; f) 

CO 

Y^Dgin-j) =D q+1 (m;t) 



Explicit form of D q (n, t) for n, t > 0: 



D (n, t) = 



(0 t<n 

Q^-p)'-"p n 0<«<f 

n < 



(Bl) 
(B2) 

(B3) 
(B4) 

(B5) 
(B6) 
(B7) 



(B8) 
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D q> o> 



Dq<a 




(BIO) 



(B9) 



D q (n, n) = D (n, n) = p' 



n 



(Bll) 
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